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^\ ' Abstract. Much use has been made of the techniques of supersymmetric quantum 

^\ , mechanics (SUSY QM) for studying bound-state problems characterized by a super- 

potential 'p{x). Under the analytic continuation ip(x) i^p{x), a pair of superpartner 
^ bound-state problems is transformed into a two-state level-crossing problem in the 

^ , continuum. The description of matter-enhanced neutrino flavor oscillations involves a 

' level-crossing problem. We treat this with the techniques of supersymmetric quantum 

mechanics. For the benefit of those not familiar with neutrino oscillations and their de- 
scription, enough details are given to make the rest of the paper understandable. Many 
other level-crossing problems in physics are of exactly the same form. Particular atten- 
tion is given to the fact that different semiclassical techniques yield different results. 
^ . The best result is obtained with a uniform approximation that explicitly recognizes the 

' supersymmetric nature of the system. 

I \ 

^ ! 

^ ■ 1 SUSY QM and the Bound-State Problem 

o ■ 

. 

\ Starting with a superpotential f{x), one can generate two superpartner poten- 

■ tials 

y± = ^2(a:)±-L^'(a;). (1) 

. .... 

O . For these two potentials the two corresponding Schrodinger equations are 



X 



2m o x 



W±{x) = E^±{x) . (2) 



It can be shown that the eigenspectrum of the system can be obtained by 
shifting the quantum numbers n of the "— " system by n — > n— 1, with the ground 
state of the "— " system discarded. That is, the spectra of the two systems are 
identical except for a single state. In applications, this property is exploited in 
the following way. Given a potential V{x), one attempts to find a superpotential 
(p{x) that will generate V{x) via Eq. (1), with one or the other sign. If this can 
be done, one can immediately generate the superpartner potential by choosing 
the opposite sign in Eq. (1). In some fortunate circumstances, the equations of 
motion for the second system are much easier to solve than the first. See Schwabl 
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(1995) for more introductory material, and Cooper et al. (1995) and references 
therein for active areas of research. 

The application of supersymmetric quantum mechanics to the solution of 
bound-state problems has been extensively developed. There has been particu- 
lar interest in semiclassical techniques. A direct primitive semiclassical (WKB) 
approach to the Schrodinger equation yields the usual Bohr-Sommerfeld quan- 
tization condition: 



where Xi and X2 are the turning points (zeros of the integrand). Since the super- 
partner potentials depend explicitly on h, the Schrodinger equations in Eq. (2) 
have a different dependence on ti than the usual case. A primitive semiclassical 
(WKB) solution of Eq. (2), which however explicitly recognizes the supersym- 
metric nature of the system [Comtet et al. (1985)], therefore yields a modified 
quantization condition: 



where xi and X2 are the turning points (zeros of the integrand). This modified 
quantization condition is exact for many systems (see Hruska et al. (1997) for a 
catalog of results), but not all [DeLaney and Nieto (1990)]. The WKB wave func- 
tions arc singular at each turning point. It is possible to avoid this problem by 
using a uniform approximation - one that is valid for all x, including at the turn- 
ing points. The approach is similar to WKB, but there is a special construction 
to cancel the turning-point singularity. A uniform semiclassical solution [Fricke 
et al. (1988)] of Eq. (2) recovers the modified quantization condition. However, 
it gives a much better wave function, which will be essential for the transition 
probability derived below. 

2 Introduction to Neutrino Oscillations 
2.1 Vacuum Neutrino Oscillations 

It is possible that the flavor and mass cigcnstates of neutrinos are not identical. 
Throughout this paper, we consider mixing between only two flavors, electron 
and muon. Then a general state can be written in the flavor basis: 




(3) 




(4) 



Ht)) 



(5) 



or the mass basis: 




(6) 



i^fi) = amplitude for the neutrino to have flavor e (/x) 
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^1 {^2) = amplitude for the neutrino to have mass mi (1712) ■ 

The amphtudcs arc taken to be time-dependent, and the kets to be time- independent. 
For two flavors, the flavor and mass bases must be related by a simple rotation. 
This rotation is taken to be between the amplitudes, with the kets held fixed, 
and is given by: 



(7) 







cos 6y — sin 9y 










sin 6.U cos Oy 







where 9y is the vacuum mixing angle. 

In the mass basis, the Schrodinger equation is 



where 
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. -^2. 





(8) 



(9) 



and is diagonal by definition. Since the neutrino masses mi and TO2 are presumed 
small, we make an ultrarelativistic expansion (using c = 1 units here and below). 
Defining E to be the common energy, and defining the mass-squared difference 
as 

dmn^ = m| — , 
one can show that -ffmass is given by 



(10) 



E- 



ml + m2 
4B 



"1 


0" 


6m^ 




0" 
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^ AE 





1 



(11) 



Below, the term proportional to the identity matrix will be dropped. 

Now we change to the flavor basis, using the rotation matrix above that 
relates the two bases. In the flavor basis, the Schrodinger equation is 



ih 



d_ 
di 







>e(t)" 




-^flavor 











where 



-fffla 



5m 
'4E 



cos 26y sin 29^ 



(12) 



(13) 



sin 29y cos 29y 

Using the rotation matrix, the amplitude to be of the electron type at a time t 
is 

^eit) = cos9y<I'i(t) + sm9y^2{t) ■ (14) 



4 J. F. Beacom and A. B. Balantekin 



The time evolution of the mass amphtudes is trivial, so that just a phase relates 
the amplitudes at a point t to those at the point t = 0. The mass amplitudes 
at t = can be expressed in terms of the flavor amplitudes at t = by use 
of the rotation matrix. Taking as initial conditions ^'e(O) = 1, = (an 

electron- type neutrino produced at i = 0), 

^ , , f ^"^"^ t\ . 2 ^ / ^fn'^ t\ 

We{t) = COS 6'„exp I +*^^^ ) +sin c'^exp I ) • (15) 

The probability for the neutrino to be of the electron type can then be shown 
to be 

P{Ue Ue){t) = \^e{t)f = 1 - sin^ 26^ sin^ (TTi/iosc) , (16) 

where the oscillation length is Lose = 4:'jTEh/Sm?, so called since this is the 

separation between extrema in the survival probability. Since the initial flavor 
state was not a stationary state of the Hamiltonian, the probability to be of 
either flavor oscillates. Note that if either the vacuum mixing angle or the mass- 
squared splitting is small, the effects of the oscillations are minimal. If either the 
source or detector has a finite size of order Lose or greater, then after averaging 
over that region, 

{P{Ue ^ Z/e))avg = 1 " ^ si^^ '^^v , (17) 

which is independent of energy. 



2.2 Matter-Enhanced Neutrino Oscillations 



In this section we assume the vacuum oscillations discussed above. Here we 
consider that the neutrino is traveling through a medium with a varying density 
of electrons, e.g., the sun. We consider the scattering 



(18) 



where Vx indicates cither flavor. In such reactions, a neutrino of a given flavor 
enters and a neutrino of the same flavor leaves. Electron neutrinos undergo both 
charged- and neutral- current reactions with the electrons, but muon neutrinos 
undergo only neutral-current scattering. For elastic forward scattering, the co- 
herence of the neutrino "beam" can be maintained. The effect of the medium is 
to modify the dispersion relation (refractive index) of the neutrino. Equivalently, 
one can say that the neutrino masses are modifled in the medium, and that the 
modification is different for electron and muon neutrinos. Below, we neglect a 
term proportional to the identity, and consider only the difference between the 
electron and muon neutrinos. Then the Hamiltonian in Eq.(13) becomes 



-^flavor — 



5w? 



C(t)-cos26'^ sin26>.„ 
sin 26*^ -Q{t)+cos29y 



(19) 
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where ({t) is related to the electron number density Ne{t) by 

2^/2GFiVe(^) 



at) 



5'm?/E 



(20) 



and Gp is the Fermi coupling constant. 

When C,{t) 0, the vacuum mixing case is recovered. When ({t) — > oo, the 
Hamiltonian is also considerably simplified. However, when ^(t) = cos2^^^,, the 
effect of the off-diagonal coupling is maximal. This is known as the resonance 
point. In a medium with a varying density, such as the sun, it is possible for an 
electron neutrino to be produced at a density above the resonance density. Then 
during its passage out of the sun it will pass through the resonance density. That 
has profound consequences for the probability that it emerges in vacuum as an 
electron neutrino. 

The flavor-basis Hamiltonian can be instantaneously diagonalized with a ro- 
tation matrix, the angle of which is called the matter angle 6{t). The bases are 
related as: 



(21) 



As the density varies, so does the matter angle. At high densities, the matter 
angle 6 7r/2; at resonance, 9 = 7r/4; and at low densities, 6 9y, the vacuum 
mixing angle (this is explained in Balantekin and Beacom (1996)). 

The instantaneous eigenvalues of the flavor-basis Hamiltonian are (reintro- 
ducing the term proportional to the identity matrix): 







"cos6»(i) 


-sin6»(t)" 










sin0{t) 


COS0{t) 







ml +1712 + Sm'^C{t) 5m 
E H T 



AE 



sm^ 20„ 



+ (C(t)-cos26'„)2. (22) 



If the density changes slowly (adiabatically), then the instantaneous diagonal- 
ization will be almost exact and the mass eigenstates will be almost stationary 
states. (Note that the Hamiltonian is exactly diagonal in the mass basis for any 
constant density.) If the density changes quickly (nonadiabatically) , then the 
neutrino can "hop" from one eigenstate to another in the neighborhood of the 
resonance (where the splitting between mass eigenvalues is at a minimum). The 
probability of this occurring is called the hopping probability Phop- This occurs 
proportionally to the extent that the instantaneous diagonalization fails. The 
instantaneous eigenvalues of -ffflavor as a function of density are shown schemat- 
ically in Fig. 1. 

Consider an electron neutrino created at a high density and a time t = 0. 
Since the matter angle will be near 7r/2, then !Z'e(0) « !?'2(0). If the density 
changes slowly, the neutrino will remain in this eigenstate. Then, if the vacuum 
mixing angle is fairly small, then 'F2{t) ~ !^/j(^) in vacuum. In this case, the 
electron neutrino created in the solar core will emerge as a muon neutrino, and 
it will be missed by experiments counting electron neutrinos from the sun. On 
the other hand, if the density changes rapidly the neutrino may emerge in the 
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Instantaneous eigenvalues of 



E + m^/2E 



E + mj/2E 



upper curve: + eigenvalue 
lower curve: - eigenvalue 

adiabatic ^^^^>^ 


resonan 

■* low density 


ce density 

high density > 



electron number density 



Fig. 1. A schematic illustration of how the instantaneous energy eigenvalues of -ffflavor 
vary as a function of electron number density. The details depend upon the mixing 
parameters, the energy, and the density profile. Beyond the left edge of the figure, 
the density is assumed to vanish; the eigenvalues are constant in vacuum and have 
the values indicated. In the figure, the initial state is assumed to be a mass = m2 
eigenstate. The final state may be the same eigenstate (adiabatic case) or the other 
eigenstate (nonadiabatic case). In the figure, the initial and final densities appear to 
be equally far from the resonance; that is not true in general. 



other instantaneous eigenstate. Then, if the vacuum mixing angle is fairly small, 
then ~ !f'e(i) in vacuum. In this case, the neutrino emerges as an electron 
neutrino after all. 

In a later section, this problem is reformulated so that these extreme limits 
of the matter and vacuum angle are not necessary. In general then, whether the 
neutrino emerges from the sun as the electron or muon type is a complicated 
function of the mixing parameters, the energy, and the form of the density pro- 
file. The probability of it emerging as the original flavor is called the survival 
probability P(fe — > fe); the general form is given in Eq. (31). 

For further introductory reading on neutrino oscillations, see Boehm and 
Vogel (1992). For a recent review of neutrino astrophysics (including oscillations), 
see Balantekin (1997). 
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3 Supersymmetric Chciracter of the Level-Crossing 
Problem 

The Hamiltonian of Eq. (19) is a typical form for a level-crossing problem. Away 
from the resonance point, the off-diagonal elements may be neglected. The di- 
agonal elements approximate the eigenvalues, and are allowed to vary. If the 
off-diagonal elements were exactly zero, the eigenvalue trajectories would cross 
(become equal) at the resonance point. This crossing will be avoided if the off- 
diagonal elements have any nonzero value, no matter how small. In the most 
general level-crossing problem, the off-diagonal elements woiild also bo allowed 
to vary. Here they are taken to be constant. Nevertheless, that is not a large 
practical restriction. Approximating the off-diagonal elements as constant is rea- 
sonable as the resonance region (the only region in which they contribute signifi- 
cantly) is usually very narrow. The explicit representation of the supersymmetry 
in Eq. (19) is given in Balantekin et al. (1988). 

Before proceeding further, we switch to working with dimensionless quanti- 
ties. We define a length scale 



L = 



hx 



5m? /AE 



(23) 



and use this to define x = t/L. In Section 1, x was used to denote a generic 
coordinate (with dimensions of length). Here and below, x is dimensionless. 
Since we will be making a semiclassical expansion, we need to be able to keep 
track of formal powers of h. For each h in the problem, we write A and consider 
A to be formally small; this is equivalent to saying that the length L is small. 
We will make expansions in powers of A, truncating the higher orders. At the 
end of the calculation, we will set A = 1. For notational convenience, we redefine 
the flavor-basis Hamiltonian as follows: 



^^dx 



avor 



where 

We have defined 
and 



■H'flavor(a^) = 



r]ip{x) \fA 
\fA —r]ip{x) 
r](p{x) = C{x) — cos 20^ 



(24) 



(25) 



(26) 



A = sin^ 26^ . (27) 

(The definition of A was misprinted in Balantekin and Beacom (1996)). The 
scaled electron density is 



ax) = 



2y/2 GpENeix) 
5w? 



(28) 
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Note that there are notation changes from previous related works [Fricke et aL 
(1988), Balantekin et al. (1988), Balantekin and Seger (1991)]; here we have 
made A and ip dimensionless. The factor r] (taken to be ±1), is introduced above 
to control the analytic behavior of the function ip{x) in the complex plane, as 
explained in Balantekin and Beacom (1996). In the expressions with ip"^ below, 
we drop r;^ = 1. 

The coupled first-order equations of Eq. (25) can be decoupled to yield 

-A^^^gM + ^\x) + iXrj^'ix)] Mx) = (29) 

and 

-A^^^JI^ -[A + cp\x) - iXw'ix)] %{x) = , (30) 

where ip{x) and A are defined in Eqs. (26) and (27). Such a simple decoupling 
is not possible in the mass basis. Eqs. (29) and (30) are explicitly of the super- 
symmetric form. This form follows directly from the (rather general) form of the 
level-crossing Hamiltonian above. The two levels are superpartners. 

These Schrodinger-like equations are similar to those for non-relativistic par- 
ticles in the presence of a complex barrier, and for convenience we use the lan- 
guage of wave mechanics to describe them. In particular, to the extent that we 
can ignore the imaginary terms in the potential, these correspond to particles 
above a barrier (since yl > 0). There are two caveats regarding discussing this 
as a barrier penetration problem. First, that our boundary conditions do not 
correspond to the usual picture of incident, reflected, and transmitted waves; in 
general, there are waves moving in each direction on each side of the barrier. 
The boundary conditions make it a level-crossing problem instead of a barrier 
problem. (The condition of a pure electron neutrino at the initial point requires 
!Z'e = 1, 'J'/i = there.) Second, the pure imaginary terms in the potentials play 
an extremely important role here, even in the asymptotic regions. These terms 
are needed not only to represent nonadiabatic transitions, but also to allow the 
local matter angle to change. Since our problem is a level-crossing problem, 
the quantity of interest is not a reflection or transmission coefficient, but rather 
P(i'e — * i^e) = \^e{x oo)Pi the probability of the neutrino being of the electron 
type far from the source. 

4 Semiclassical Solution of the Equations of Motion 

4.1 General Form of the Solution 

In general, the survival probability has terms which depend on the source and 
detector positions. In this paper, these interference terms are considered to be 
averaged away by the finite source and detector sizes (or by a varying distance 
between them, such as that due to the motion of the earth) . The remaining term 
below must also be considered to have been appropriately averaged. 
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The general form of the survival probability is 

P(z.e ^ /^e) = i [1 + (1 - 2Phop) COS 29^ COS 20^] . (31) 

The two matter angle terms account for the rotations between the mass and 
flavor bases at the initial and final (in vacuum) points. (The matter angle thus 
ranges from tt/2 at infinite density to Oy in vacuum. At the resonance, 6 = 7r/4.) 
Phop is the probability of hopping from one mass eigenstate to the other during 
the transit through the resonance region. This result for the survival probability 
is completely general. The heart of the problem is to determine Phop- Various 
semiclassical results for Phop are discussed below. The variation of the averaged 
survival probability with Sm^/E for two choices of the vacuum mixing angle is 
shown in Fig. 2. 

The general shape of the survival probability vs 6m^/E is that of a pit. High- 
and low-energy electron neutrinos have a relatively high probability to remain so, 
whereas medium-energy electron neutrinos have a relatively low probability to 
remain so (and hence a relatively high probability to transform into muon neutri- 
nos). The solar neutrino detectors are primarily sensitive to electron neutrinos. 
The energy-dependent suppression of electron neutrinos as in Fig. 2 explains the 
solar neutrino data well. 

4.2 Primitive Semiclassical Solution 

A primitive semiclassical (WKB) solution can be made to Eqs. (29) and (30). 
It can be shown that this is exactly equivalent to the adiabatic solution. The 
adiabatic solution is obtained when the Hamiltonian changes so slowly that the 
instantaneous diagonalization always holds with high accuracy. That is, there 
are no transitions between mass eigenstates, so Phop = by definition. This 
result is not shown in Fig. 2. If it were, it would match the exact result in all 
but the left-hand rise of the pit. Where the other curves rise up, the WKB result 
would remain flat, at the level of the bottom of the pit. 

In the bound-state problem, even the primitive semiclassical solution gave 
an excellent answer for the quantization condition. As noted, in the neutrino 
problem this approximation does not predict Phop- The problem comes from the 
well-known singularity of the WKB wave functions near the turning points. In 
this problem, the turning points are in the complex plane, near the resonance 
point (which is on the real axis). In the bound-state problem, the turning point 
singularity is not as crucial since the quantization condition only counts nodes 
between the turning points. In this problem, we need the wave function itself to 
be accurate, since its modulus squared determines the survival probability. 

4.3 Landau-Zener Solution 

Far from the resonance, the propagation is adiabatic, and can be trivially solved 
for any density profile. A reasonable approach to simplifying the problem is to 
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Average survival probability 




Sin'ie^ = 0.7 

0.0 _in _q _o _7 _£ _c A _T 

10 10 10 10 10 10 10 10 
dia/E (eV'/MeV) 

Fig. 2. The electron neutrino survival probability vs. the mass-squared difference par- 
rameter for two different vacuum mixing angles. The solid line is given by the method 
of Balantekin and Beacom (1996). The dashed line is the exact (numerical) result. The 

dotted line is the hnoar Landau-Zcnor result. In the top figure, the hues are indistin- 
guishable. An exponential density with parameters chosen to approximate the sun was 
used [Bahcall (1989)]. The region leftward of the lower left corner of the trough is the 
nonadiabatic region. 

approximate the density profile as linear in the resonance region. This is the basis 
of the Landau-Zener solution. In the flavor-basis Hamiltonian of Eq. (25), this 
gives the usual Landau-Zener setup: linear variation of the diagonal elements, 
and constant off-diagonal coupling. In order to make Fig. 1 more general, the 
variation is shown versus the density itself (instead of x). 

The solution is obtained from the decoupled equations (29) and (30). With 
({x) X, the differential equations can be reduced to the defining equation 
for the parabolic cylinder (Weber) functions. Given the initial conditions = 
IjiJ'p = 0, the solution for ^e{x) at any point is straightforward. (Actually, 
the boundary conditions require some care since the linear density eventually 
becomes negative, which is unphysical.) Prom the asymptotic form appropriate 
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far after the resonance, one can extract the hopping probabihty Phop from the 
the survival probabihty P{Pe ^e) = We{x — + +oo)p. The result is 



Phop = exp 



(5m^ sin^ 29.u 
UeK cos26I„ 



at) 



(32) 



Because the Landau-Zener solution is based on the exact solution for the 
linear density, there are no turning-point singularities as with the WKB solu- 
tion. Thus the approximate wave function and hence Phop ai"© fairly reasonable. 
However, there is some inaccuracy due to the fact that a general density does 
vary more than linearly through the resonance region. The fact that the Landau- 
Zener result has the right general behavior but is not very accurate can be seen 
in Fig. 2. 



4.4 Uniform Semiclassical Solution 

In this section we outline a uniform semiclassical solution to the problem [Bal- 

antekin and Beacom (1996)]. An arbitrary monotonic density profile is allowed, 
as are nearly arbitrary mixing parameters. We start with Eq. (29): 

-A^^^^M _[A + cp^x) + iXvif'ix)] Mx) = . (33) 

There is no known exact solution of this for arbitrary C,{x). Recall that A is a 
formally small perturbation parameter (we set A = 1 later), and that r/ = ±1. 
Compare this to 

-A^ -[n + S'± iXv] UiS) = , (34) 

which is solvable in terms of parabolic cylinder functions. Identification of the 
turning points of each of these will be crucial. At lowest order in A, they are the 
points for which A + ip'^{x) =0 and f2 + S'^ = 0, respectively. The turning points 
are complex conjugate pairs. 

Suppose wo could find a change of variables S ~ S{x). Using this, we could 
deform the simple parabolic "barrier" of the second equation into the more 
complicated "barrier" shape of the first. Then the solutions of the general case 
could be expressed in terms of the solutions of the simple case. Such a formal 
solution can be written as 

Mx) = -j^UiSix)) . (35) 

Starting with this formal solution, and the two Schrodinger equations above 
for ^e{x) and U{S), one can determine a differential equation for S = S{x). Un- 
surprisingly, that is a nonlinear equation, and at least as difficult as the original 
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Schrodinger equation. However, this equation can be profitably subjected to a 
semiclassical solution. We expand as 



S{x)=So{x) + XSi{x)... , 



(36) 



and truncate at second order. Separating the various orders in A determines a 
series of equations. The solution of the 0{X) equation determines So{x): 



and also /?: 



i-So{x) I px 

/ dSo^n + S^= / dt^A + ^^x) 
J7 = — / dx^/A + if^Jx} , 

Jxo 



(37) 



(38) 



where xo and Xq are the turning points (zeros of the integrand). The 0{X^) 
equation determines Si{x): 



Siix) 



IT) 



2^n + s'^{x) 



In 



In 



f2 



So{x) + ^f2 + Si{x) 



(39) 



The change of variables S = S{x) is actually a mapping of the complex 
x-plane to the complex 5-plane. In order to avoid spurious branch cut discrep- 
ancies, the mapping must be chosen to not fold or flip the plane. This is ac- 
complished with some appropriate choices of signs. Further, the turning points 
in the x-plane must be mapped onto the turning points in the 5-plane. That 
demand fixes the definition of f2 to be the one given above. 

The approximate (but uniformly valid) solution is 



f2 + S^{x) 
A + ip^{x) 



1/4 



U{So{x)+XSi{x)) 



(40) 



At the turning points, the denominator vanishes, just as in the WKB solution. 
However, because of the matching of turning points in the mapping, the nu- 
merator vanishes at the same points. That cancels the singularity and gives an 
excellent approximation to the wave function at all points, including near the 
turning points (and hence near the resonance, which is on the real axis near the 
turning points). 

As noted above, the functions U{S) are parabolic cylinder functions. Using 
the defining equations for So{x) and Si{x) given above, one can solve for i^e{x) 
for any x. The implicit definition of 6*0(2;) makes it analytically solvable only 
for large |x| and large \So\, which holds far from the resonance. However, both 
the initial and final point can be taken to be far from the resonance. At each 
point, the general solution can be written in terms of two independent parabolic 



A Semiclassical Approach to Level Crossing in SUSY QM 13 

cylinder functions with arbitrary coefficients. After taking the asymptotic forms 
(meaning well away from the resonance), one has: 

^e{x -oo) (41) 
= C-Cos9{x)ex.p{+iIp{x,Xi)/X) + D-sm9{x)ex.p{—iIp{x,Xi)/X) , 

^e{x +00) (42) 
= C+sin^(a;)exp(— i/p(x, a;i)/A) + Z)_|_ cos^(a;) exp (+i/p(a;, a;i)/A) 
= [ci sin^i exp (+iRe7p(a;i, a;o)/A) + C2 cos^, exp (— iRe7p(a;i,a;o)/A)] 

X sm9{x) exp (— iRc/p(.x, X())/X) 
+ [cl cosOi exp {—iRcIp{xi, xq)/X) — sinOi exp {+iReIp{xi, Xo) / X)] 
X cos0{x)exp{+iReIp{x,xo)/X) , 



where 



Ip{x,Xi)= I dx^A + ip^{x). (43) 

J Xi 



Above, x is a general point, Xi is the initial point, and x^ is a tinning point. Now 
one must solve the problem of connecting the coefficients C_ and D- (known 
from the initial conditions) to the coefficients C+ and 1?+ (or equivalently Ci 
and C2). Tedious algebra reveals the solution to the connection problem to be: 

c. = .£^n-n,.,r„. (^^\^^,^^ 2.sin(..) (44) 



V 2 y ^/2 " V 2 
C2 = e-'"'' . (45) 

Comparison of \^e{x ^ +00) I ^ above to the general form 

P{ve ^ye) = \[l + {l- 2Phop) cos 29, cos 29^] (46) 

reveals that: 

l-ftop = |ci|2 = l-e-^", (47) 
P^.op = \c2? = e-"\ (48) 

where 

2i f 5rn?\ I 

j dt^&m^ 29, + (C(t) - cos29,f . (49) 

The dimensions have been restored and A and rc'placed by their definitions; re- 
call that C(t) ~ density and Iq, are the turning points (zeros of the integrand). 
For a linear density, the Landau-Zener result is recovered, i.e.. 



A r, sin^ 29, 

linear density : il = —— — 

^ AEn cos29. 



at) 



(50) 
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However, any arbitrary monotonic density may be used. For example, for an 
exponential density, 

-1 

(51) 

which is the leading term in the exact result. That the approximation is rather 
good over a wide variety of parameters can be seen in Fig. 2. 

5 Conclusions 

In a semiclassical approximation, one always finds that the wave functions have 
terms like 

!Z^~exp(( )/n) , (52) 

where ( ) indicates a phase integral and various numerical factors. There is an 
essential singularity in tf' in the formal limit /i — > 0. That is, ^ does not have 
a well-defined value in this limit. Difi'erent approximations (i.e., different ways 
of limiting ?i ^ 0) therefore give different results. Not all semiclassical approx- 
imations are equivalent; other criteria must be used to decide which approach 
to use. In this level- crossing problem, the best results come from a uniform 
approximation (valid for all x, including the turning points) that explicitly rec- 
ognizes the supersymmetric nature of the potential. Given the success of similar 
approximations in the bound-state problem, this is perhaps to be expected. 

Even the uniform supersymmetric approximation eventually breaks down. 
In the original problem (in the figure above, based on the exponential density), 
we considered only the two turning points closest to the real axis. In general, 
there are more, further out in the complex plane. The distance of the the turning 
points from the real axis scales like |C/CI (and for the exponential density, integer 
multiples of this). The analog potential (based on the linear density), only has 
two turning points. This mismatch in turning-point topologies eventually makes 
the mapping function multivalued, at which point the approximation breaks 
down. The problem is worst when the length scale Eh,/5m^ becomes large and 
"sees" further into the complex plane. 

This is an interesting practical problem, in part because there is good ev- 
idence for matter-enhanced neutrino oscillations. The technique given here is 
applicable for any monotonic density profile (and not just in the sun), and for a 
large range of mixing parameters [Balantekin and Beacom (1996)]. In addition, 
the techniques here should also be useful in a variety of the continuum level- 
crossing problems that are ubiquitous in quantum mechanics. As noted, certain 
kinds of two-level systems arc automatically of the supersymmetric form. The 
techniques of semiclassical SUSY QM for the bound-state problem can be carried 
over, as shown in our example. 



5m? t{t) 
exponential density : i? = (1 — cos 20^) -— - 

AhjJi Q\t) 
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